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Abstract 

Q ■ We consider front propagation in a reactive Boussinesq system in an infinite vertical 

strip. We establish nonlinear stability of planar fronts for narrow domains when the 
Rayleigh number is not too large. Planar fronts are shown to be linearly unstable with 
respect to long wavelength perturbations if the Rayleigh number is sufficiently large. 
We also prove uniform bounds on the bulk burning rate and the Nusselt number in the 
KPP reaction case. 



> . 

^ ; 1 Introduction 

o : 

The presence of fluid flow may have a strong effect on reaction processes in many situations 



> 



! of interest |L9|, |28| . Many recent studies focused on the influence of passive advection on 

combustion. During passive advection the reacting material is carried by a flow, but the 
flow is not influenced by the reaction. Traveling waves and pulsating fronts in periodic flows 



have been studied in || |], ||, |25], |26| ; the thin front limit was considered in ]TT|, [12], |17]] using 



homogenization techniques. Estimates on the speed of front propagation for different classes 
of flows have been derived in [[l], |5|, [16], || [13[] . Numerical studies of the effect of a passive flow 
on combustion were carried out in |TJ|, ^, |S], |2"3|. This list is far from complete; we refer 
to the recent reviews [[[], [27] for further references. 

However, there have been few rigorous works on models with feedback of the chemical 
reaction on the flow field. In this paper we study a simplified active combustion model 
[18] , |22| , gy , |2l| in which the reaction does influence the flow. The feedback of the flame 



on the fluid is taken in a Boussinesq approximation. The model couples thus an advection- 
diffusion-reaction equation for the temperature with an incompressible Navier-Stokes system 
driven by temperature differences. We study this problem in a two dimensional strip of 
infinite vertical height and finite horizontal width. The vertical direction is the direction of 
gravity. The system admits planar fronts as particular solutions. These fronts correspond 
to traveling solutions of the one dimensional reaction-diffusion system without horizontal 
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variation. We study them in the context of the larger reactive Boussinesq system. Coupling 
with the Boussinesq system introduces at least two new interesting effects. One is symmetry 
breaking: gravity breaks the vertical symmetry of the reaction diffusion systems. This has a 
dynamical effect: fronts connecting low regions of hot fluid to high regions of cold fluid are 
susceptible to the Rayleigh- Taylor instability [|7| . The second effect is due to the introduction 
of new horizontal degrees of freedom: new length scales are introduced. When one ignores the 
fluid advection, the planar fronts have a characteristic thickness 5 (|2.11|) , which is determined 



by the thermal diffusivity and the characteristic reaction time t c ( |2.13| ). Using these as length 



and time units, three significant nondimensional parameters emerge. One is the aspect ratio 
A ( 2.23 ), the ratio of the horizontal width of the strip to the thickness of the planar front. The 



second parameter is the Prandtl number a ( |2.f 9| ), the ratio of kinematic viscosity to thermal 



diffusivity. The third important parameter is the Rayleigh number p ( p.20| ) which measures 
the relative strength of the buoyancy force on the scale of the front thickness. 

We prove three kinds of results: stability, instability and uniform bounds. Stability is 
analyzed in Section [|. If the aspect ratio A is small, then the only traveling solutions are 
planar fronts, and all solutions become eventually planar. This is a consequence of the fact 
that diffusion acts rapidly across a narrow strip. This stability mechanism is quite robust and 
operates for all kinds of nonlinearities. If the Rayleigh number is large enough, then the planar 
front loses stability to long wave perturbations, which are present if the aspect ratio is large 
enough. The instability is of Rayleigh- Taylor type and is also expected to be quite robust. We 
note that nonplanar traveling waves have been recently shown to exist (|2(| ^TJ) in reactive 
Boussinesq systems with bistable nonlinearity at large enough Rayleigh numbers. Our results 
on planar front instability agree qualitatively with recent numerical results [22]|. The proof 



appears in Section and is based on the intuitive idea that, in the presence of gravity, the 
neutral mode - corresponding to the broken vertical translation symmetry - misaligns over 
long horizontal distances, and gives rise to a long wave unstable mode. Technically, one needs 
to deal with a variable coefficient linear operator in which the instability of heavy fluid on top 
of light fluid is exploited using the monotonicity of the front. In Section |3] we prove bounds 
for arbitrary solutions of front type. The results state that the bulk burning rate, speed of 
front, gradients of temperature, and fluid quantities, are all bounded. The bounds do depend 
on the aspect ratio. The proof applies only to concave KPP nonlinearities. The results imply 
that there are no accelerating fronts in this system. The proof of the upper bound passes via a 
lower bound: the temperature gradients squared bound from below the bulk burning rate 0. 
On the other hand, the bulk burning rate is bounded above by the sum of laminar burning 
rate and temperature gradients to power one. This implies that the temperature gradients 
are bounded, and the rest follows. 

2 Reactive Boussinesq fronts 

The reactive Boussinesq equations are 

dv 

— + v ■ Vv + Vp - uAv = gATe z , (2.1) 
V-w = 0, (2.2) 
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8T v 2 

— +v-VT-nVT=^f{T). (2.3) 
at 4k 

Here v (x, t) = (u(x, t), w(x, £)) and T(x, t) are the velocity and the (normalized) temperature. 

The vector e z is the unit vector pointing in the direction opposite to the pull of gravity. The 

constants v, k, A, Vq, g are all positive, and represent, respectively, the kinematic viscosity, 

thermal diffusivity, thermal expansion coefficient multiplied by the temperature variation 

scale, speed of purely reactive-diffusive front and acceleration of gravity. The variables x = 

(x, z) belong to a strip, with z6K and x G [0, L\. The boundary conditions for the normalized 

temperature are front conditions, 

T(x,z,t)—*\ as z — > — oo, T(x, z, t) — > as z — ► oo. (2.4) 

We assume throughout the paper that the initial condition T (x, z) satisfies < T (x, z) < 1. 
The velocity vanishes at the two ends of the strip: 

v(x, z, t) — > as \z\ — > oo (2.5) 

The lateral boundary conditions are periodic: 

T(x + L, z, t) — T(x, z, t), v(x + L, z,t) — v(x, z,t). (2.6) 

We consider the vorticity 

dw du . 

where u is the component of v in the x direction and w is the component of v in the z direction. 
The momentum equation ( |2.1|) implies 

duo dT 

— + v-Vu- uAu = gA— (2.8) 
at ox 

The nonlinearity f(T) satisfies /(0) = /(l) = 0, is a smooth function and is assumed to be of 
one of the following three types: bistable, ignition type or concave KPP. Bistable nonlinearities 
are such that there exists 9 G (0, 1) so that f(T) < for T G (0, 9) and f(T) > for T G (9, 1). 
Ignition type nonlinearities are such that there exists 9 G (0, 1) so that f(T) = for T G [0, 9} 
and f(T) > for T G (9, 1). Concave KPP type nonlinearities have f(T) > for T G (0, 1) 
and f"(T) < —M for T G [0, 1] with M > 0. In their case we assume that the slope at zero 
is positive and normalized, /'(0) = 1, while the slope at 1 is negative, f'(l) = —m < 0. We 
will use the KPP properties only in the proof of a uniform upper bound for the burning rate 
presented in Section || We summarize the essential properties of the KPP reaction used in 
that proof: 

< f(T) < T, -f"(T) > M > 0, —f{T) < m, for < T < 1. (2.9) 
The Boussinesq system has flat traveling wave solutions 

T fr = r{z -ct),v fr = 0. (2.10) 



The momentum equation ( |2.1|) holds in this case because the pressure can balance a temper- 
ature that depends on z and t alone. The speed c takes all values c > v o in the KPP case, 
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and is unique in the bistable and ignition case. The profile t(z) is monotonically decreasing 
in all three cases, and obeys 

kt" + ct> + ||/(r) = 0, 

dr 

where r' — — . We choose now dimensional units. First we note that the function r varies 

dz 

on length scale of the order 

*=-, (2-11) 

that is, the function r has the form t(z) = P(z/S), where P obeys the equation 

P" + 2cP' + f(P) = 0. (2.12) 

Here c = c/vq is a nondimensional constant that is determined solely by the form of the 
nonlinearity f{T), and that in turn determines the traveling front speed c = cvq in the 
original dimensional variables. We choose reaction units: we take fl2.11|) for length unit and 

we take 

Ak 

f e = -j (2-13) 
for time unit. The velocity fluctuations are written as 

«<*■«> =Ht" 9 (2 - l4) 

and the temperature fluctuations are written as 

6(x 1 t) = §(j^y (2.15) 

Using these units, rescaling, using x = (x,z) = (x new ,z new ) = (x oid /5, z old /5) and t = t new = 
toid/tc, and dropping tildes, we derive the nonlinear equations 

— + v ■ Vu - aAu = ap— (2.16) 
at ox 

and 

^ + v-VT-AT = f(T), (2.17) 

where v = (u, w), with 

dui . dui , 

Al, = -&' A »=& < 2 - 18 > 

which follows from (|2.7| ) and incompressibility. The nondimensional parameters are the 
Prandtl number 

a = - (2.19) 
and the Rayleigh number across a laminar front width 

p = - . 2.20 

KV 
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The boundary conditions in z are of front type 



T(x,z,t) - 
T(x, z, t) - 
v(x, z, t) - 
u)(x, z, t) 



1 as z — ► — oo, 
as z — > +00, 
0, as \z\ — > 00. 
* as Id — > 00. 



(2.21) 



The boundary conditions in x are periodic 

T(x+ A, z, t) — T(x, z, t), v (x + A, z, t) = v(x, z, t), lu(x + X,z,t) — ui(x,z,t) 
with period 

3 General bounds in the KPP case 



(2.22) 
(2.23) 



We study general solutions of the system (|2.16| , |2.17| , |2.18|) with a KPP nonlinearity that 
satisfies assumptions ( |2.9| ), and with front-like initial data: functions that approach the values 
as z — > 00 and, respectively 1 as z — > —00 at least at an exponential rate. We denote 
D = [0, A] x E. We will use the notation 



\9\\h 



A 



\g(x, z)\ 2 dxdz 



(3.1) 



D 



for the normalized L 2 norm. We consider the bulk burning rate 



dT(x, z, t) 



dxdz 



(3.2) 



as a measure of the reaction rate. Using equation ( 2.17 ), incompressibility of u and the 
boundary conditions for u and T as z —>■ ±00 we have also 



(3.3) 



V(t) = - / f{T(x,z,t))dxdz 
A Jd 



The time average of V is denoted V 



— 1 f 

V(t) = - V(s)ds. 

t Jo 



(3.4) 



Lemma 1 Assume that there exists a constant a G R so that the front-like initial data T (x, z) 
obeys 

Tq(x, z) < exp (a — z) (3.5) 



and 



1 - T (x, z)) < exp (a + z). 



(3.6) 
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Then the solution of obeys the bounds 



T(x, z, t) < exp 



a — z + 2t + / \\w(', s)\\ Loads 



and 



a + z + t 



(1 — T(x, z, t)) < exp 
for allt > 0. 

Proof. For the bound of T we seek a supersolution of the form: 



w(-, s) \\L°°ds 



(3.7) 



(3.f 



9 + (z, t) = exp [— az + a / \\w(-, s)\\ L °cds + bt + a], a > 0. 

Jo 



Using fl2.9Q we get 



+ A6 + -f(6 + ) >0 



if 6 > a 2 + 1. We chose for simplicity of exposition a = 1, and took the most economic 6 = 2. 
For the bound of 1 — T we seek a subsolution for T of the form 

6L(z, t) = 1 — exp [az — a / \\w(-, s)\\L<*>ds + bt + a] 



and, using the fact that / > on [0, 1], the condition 

86. 



dt 



follows if 



+ v ■ V0_ - A0_ - f(9-) <0 



b>a 2 . 



Again, we chose a = 1 for simplicity, and took the best 6=1. The proof shows that similar 
inequalities hold if the exponential rate of decay at infinity of the initial data is different than 
a = 1, or if the initial data approaches its limit at negative infinity at a different exponential 
rate than a — 1. 

Let us consider the average quantities 



and 



W(t) 



N(t) 



t 



\w(; s)\\ L °°ds 



||VT(-,5)|| 2 i2 ^ 



(3.9) 



(3.10) 



representing the average maximum vertical velocity and temperature gradient squared. 

Lemma 2 Consider front-like initial data that satisfy §3. dp . Then the solutions of ($■!% ) 
obey 

V(t) < W(t) + 2 + - (3.11) 

t 

for all t > with 7 depending on the initial data. 
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Proof. First we write 



V(t) 



xt 



(T(x, z, t) — T (x, z)) dxdz. 



D 



which we bound as 



V(t)< 



Xt 



dx 



(1 — T (x, z)) dz + / T(x,z,t)dz 



(3.12) 



using the fact that T(t,x,z) < 1. Now, denoting 

B 1 (t) = a + 2t+ f \\w(',a)\\ L oods, 
Jo 

we have from ( |3.7| ) that 

poo 

/ T(x, z, t)dz < 1, 

JB^t) 

while, because T < 1, we have 



(3.13) 



/ T(x,z,t)dz < Bx{t). 
Jo 



We use these bounds in (|3.12| ) and obtain ( |3.11|) . Now we bound W{t) below by N{i): 

Lemma 3 Consider front-like initial data that satisfy ( \3. b\ ). Then the solutions of ( \2. 1 lj ) 
obey 

m + 2 . . 2m + 3 T , 

N(t) < —Wit) + — + 7 (3.14) 

with m and M given in (\2.9j) and with T depending on the initial data. 



f'(T(x, z, t)) dT ^ Z,t ^ dxdz. 



Proof. We start by computing 



>4 



D 



Using ( p.!7| ) and integrating by parts, using incompressibility of u and the boundary condi- 
tions, we obtain 



dV 1 
~dt ~ X 



f'(T(x, z, t))f(T(x, z, t))dxdz 



D 



f"(T(x, z, t))|VT(x, z, t)\ 2 dxdz. 



D 



Using (|2.9| ) we deduce 



Taking a time average we get 



1 



^ + mV(f)>M||Vr(.,f)lli»- 



- (V(t) - V(0)) + mV(t) > MN(t). 



(3.15) 



We observe that 



i r r dx r 2K ) 

V(t) = - f(T(x,z,t))dxdz = — f(T(x,z,t))dz (3.16) 

A J D JO * J-oo 

+ r x dx f Bl(t) f{T{x ^ Z}t))dz+ f X <^_ r f{T (x,z,t))dz, 

JO A J-B 2 (t) JO A JBx{t) 



where B\{t) is given in ( |3.13| ) and 



B 2 (t) = a + t+ / \\w(-,s)\\ L ™ds. 
Jo 

We use the inequality f(T) < m(l — T) which follows from ( |2.9| ). Then we have 

~ x dx r BaW 



A 



f(T(x, z, t))dz < m, 



as follows from ( |3.8|) . Similarly, 

-A 



f dx r°° 

/ ^ / f(T(x,z,t))dz<l 

Jo A JB t (t) 



follows from ( |3.7| ) and ( |2.9| ). The second term on the right in ( |3.16| ) is bounded by Bi(t) + 
B 2 (t). Thus, returning to (|3.15| ), we have 



MN(t) < mV(t) + 3 + 2W(t) + -. 



In view of (|3.11| ) of the previous lemma, (|3.14|) is proven. 

The next step consists in bounding the quantity W(t) in terms of N(t), using the vorticity 
equation (|2.16|) . 

Lemma 4 There exists an absolute constant C so that, for all t > one has 

W{t) < CA 3/2 ^pXy/N^j + -L||u;o|| z2 } , (3.17) 



where Wq(x, z) is the initial data for u>(x, z, t). 

Proof. We multiply (|2.16|) by uj and integrate. After one integration by parts we obtain 

f Hx , z , t )\^ + a [ \Vu(x, z, t)\^ = a P [ u(x, z, t) ^ M) (3.18) 
2dtJ D A J D A J D dx A 

We introduce 

T(z,t) := J T(x,z,t)-£ 

and note the obvious fact that 

dT(x, z, t) _ d(T(x,z,t)-T(z,t)) 
dx dx 
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Inserting this expression into the right hand side of ( |3.18| ) and integrating by parts on the 
right side we obtain 



1 d 
2~dt 



\u>(x, z, t)\ 



,dxdz 



A 



+ a / \Vuj(x, z, t)\ 



,dxdz 



D 

-op 



D 



dcu(x, z, t) 
dx 



(T{x,z,t)-T(z,t)) 



A 

dxdz 

~1T 



D 



Using Young's inequality together with the inequality 

,2 dxdz 



\T(x, z, t) — T(z, t) 



we deduce 



D 



1 d 



A 



< A' / \VT{x,z,t)\ 



D 



dxdz 

~1T 



2dt 



a 

< - 
~ 2 



,dxdz 



-— / \u>(x,z,t)\ — her / | Vuj(x, z, t)\ 



D 



A 



D 



,dxdz 

~1T 



D 



duj(x, z, t) 



dx 



dxdz ap 2 X 2 



A 



\VT(x,z,t)\ 



2 dxdz 



D 



A 



(3.19) 



Integrating (|3.19|) in time we deduce 



r/.s / \\/uj(x,Z,s)\ 2( ^^ < p 2 X 2 N(t) + — Hwollia 

A at 



(3.20) 



'0 JD 

Let us represent the function w in terms of its Fourier series 

Akx 



w(x,z,t)= Wk(z,t)e l 



and note that, in view of incompressibility, Wq(z, t) is independent of z, and hence the bound- 
ary conditions at z ± oo imply that 

w (z,t) = 0. 
In view of Q2.18 ), the embedding inequality 



k(-,t)||L-<CA 3 / 2 ||V^(-,t)|U 2 



(3.21) 



follows. The constant C is absolute (recall that the L 2 norm is normalized by A ( |3.1| )); it can 
be computed either by using Parseval's identity or using ( |5.16| ) below. From (|3.21| ) and ( |3.20| ) 
we deduce ( ^.17| ), using the Cauchy-Schwartz inequality for integration in time. 



Theorem 1 Solutions of $2.1(\ ), ($.1T\ ), ( \2.18{ ) with front-like initial data obey 



N(t)< 



C 2 (m + 2) 2 2 5 4m + 6 K x K 2 



M 2 



-p z \° + 



M 



Vi t 



M 



U *3 



(3.22) 
(3.23) 



and 



lim sup V(t) < 2 + 

t— >oo 



C 2 (m + 2) 



p 2 \ 5 + C 



Am + 6 



-P 



A 5/2 



(3.24) 



M V M 

wt/i M, m gwen in (2.i>), C the absolute constant of US. 21 ) and with Ki, K2, K3, K4 depending 



on the initial data . 

Proof. We insert ( p,17|) into the right hand side of ( p.!4|) . We obtain 

N(t) < 



M 



2m + 3 T 



MVat 

Now we use Young's inequality 



and deduce ( ggg ). Then (|3~23|) follows from and (CT) , while follows from ( p^3|) 

and (|3Tfl 



Remarks. 1. In terms of the original dimensional variables of ( |2.1| ),( [2~2"D ,( [2.3| ) the bound 
on the wave number (|3.22 ) means that the Nusselt number on scale L 



Nu 



lim sup - / / / I vT(x, z, s)\ 2 dxdzds = A lim sup N(t) 



is bounded as 



Nu < 



C 2 (m + 2f 2 . 4m + 6 



M 2 



with i?a the Rayleigh number on scale L: 



Ra 



-Ra z + 



gAL 3 



M 



A 



VK 



(3.25) 



(3.26) 



2. Numerical evidence [g2] shows that the V^(t) is indeed an increasing function of A and p. 



4 Nonlinear stability of planar fronts in narrow do- 
mains 



In this section we consider the reaction-diffusion Boussinesq problem ( |2.16p - ( |2.18| ) in a 
narrow domain, i.e. for small aspect ratio A. The nonlinearity / is of either one of the three 
types: KPP, ignition or bistable. We keep the time and space units chosen above fl2.ll , |2. 13 
[2.14| , |2.15| ). We prove two results. The first one concerns traveling solutions of the form 



T(x, z, t) = T(x, z — ct), v(x, z, t) = v(x, z — ct). 



(4.1) 



The result states that, if the Rayleigh number p is sufficiently small, then such solutions must 
be planar fronts. Planar fronts are solutions of the form 



T(x, z, t) = P{z — ct), v(x, z, t) — 



(4.2) 



for which T does not depend on x and v vanishes. Planar fronts do exist for c > 2 in the 
KPP nonlinearity case and for a unique front speed c* in the ignition and bistable cases. 
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Theorem 2 There exist constants C\ > and C2 > such that if A < Ci ; and p < C2/X 3 , 
then the only solutions of ( $.ldj )- ft2.1$ ), (\2.21[ \2. 22j) of traveling front type are planar 

fronts of the form j\4-2j )- 

The second result in this section is about arbitrary solutions. We show that all solutions of 
the Boussinesq system in a narrow domain become eventually planar: 

Theorem 3 There exist constants G\ > and C2 > so that if A < G\ and p < C2/A 3 , then 

\\uj(.,t)\\ L 2 + \\T x (;t)\\ L2 ^0 as t-> +00. (4.3) 
Moreover, the front speed is uniformly bounded: 

limsup V(t) < 2. (4.4) 

Proof of Theorem Let T(x, z — ct), v (x, z — ct) be a traveling front solution of 
(1218)) , then 



-cT, + v ■ VT = AT + f(T) (4.5) 
— ct> 2 + v ■ Vv + Vp = aAv + <ypTe z 
—coj z + v ■ Vu = aAu + apT x . 

We multiply the second equation in (|4.5|) by v and integrate to obtain 

o"||Vf || 2 = op J T(x, z)w(x, z)dxdz. 
In this section we will use the plain L 2 norm 

ll^ll 2 = / \g(x, z)\ 2 dxdz. 
The right side may be bounded as in the proof of Lemma f| to obtain 

a||Vw|| 2 = crp / T(x, z)w(x, z)dxdz < crp\\T — T\\ \\w\\ < apX\\T x \\ \\w\\, (4.6) 



where 

dx 



T(z) = [ T(x,z,t)- 
Jo 



A 

We used above the fact that w(x, z) has mean zero in x because of incompressibility and 
boundary conditions at z — >• 00. Furthermore, it satisfies the Poincare inequality so that 
|| io 1 1 < CA||Vf ||, and the incompressibility of v implies that ||Vf || = ||co>||. Therefore Q4.6|) 
implies that 

||w|| < CpA 3 !!^!! 

and 

\\u\\<C P X 2 \\T x \\. (4.7) 
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We differentiate the first equation in ( |4.5| ) in x, multiply by T x and integrate to obtain: 

-J T x {v x -VT)dxdz = J \VT x \ 2 dxdz- J f'(T)T 2 dxdz. (4. 

We use the Poincare inequality for T x to obtain a bound for the right side: 

J \VT x \ 2 dxdz - J f'(T)T 2 dxdz > (l - CA 2 ) J \VT x \ 2 dxdz > -\\VT X \\ 2 

if A < C\ is sufficiently small, with C\ an absolute constant that may depend only on m\ 
max <T<i \f'(T)\. The left side of (^4.8| ) is bounded by 



T x (v x ■ VT)dxdz 



T{v x ■ VT x )dxdz 



— \\ v x\\ l|VTJ| < ||a;|| II VTJ 



The last two bounds imply that HVT^H < 2||co>|| for a sufficiently small A. Then the inequality 
(|4.7| ) and the Poincare inequality for T x imply that 



\\VT X \\ < Cp\ 2 \\T x \\ < CpX 3 \\VT x \\. 
Therefore there exists a constant Co > such that no x-dependent traveling front exists if 

(4.9) 



This finishes the proof of Theorem |2|. We now prove Theorem |3|. 

Proof of Theorem [3|. We multiply the evolution equation for v(x,z,t) by v = (u,w) 
and integrate to obtain 



ld_ 
2~dl 



+ a\\ Vv (•, t) || 2 = op \ T{x,z,t)w{x,z,t)dxdz. 



(4.10) 



The right side of the above equation may be bounded as in (|4.6|) . We integrate ( 4.10| ) in time, 
and use ( |4.6| ) and the Poincare inequality for w(t, x, z) to get 



< crpA 



w(;s)\\ 2 ds<Op\ / \\T x {; S )\\\\w{;s)\\d S +\\v{;0)\\ 2 (4.11) 
(I JO 

1/2 / r t \ 1/2 







\Tj; S )\\ 2 ds) { I \\W(; S)\\ 2 ds) + \\v(;0)f 







ap\ ( C [\ 2 pA 3 



J*\\T x (.,s)\\ 2 ds^ +\\v(;0)\t 



Therefore we have 



V(-, S )|| 2 d S <C> 2 A 6 /V«Mll^fe + — IK.O)!! 1 
o Jo a 



Furthermore, once again, after integrating ( (4.10 ) in time we obtain 



t rt ft (~i 

2j„/ „\ / \\rn ( „MIIL../'„MI j„ ^ /~<„2x4 / \\rr /„M|2 j„ , ° ||„./ nM|2 



£*>(■, s) || 2 (is < pA / ||T x (-,s)||||w(s)||ds < Cp 2 X A / 11^(5)11^5+ -||w(-,0)|| 2 . (4.12) 



(7 
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Furthermore, we differentiate the first equation in ( 2.17| ) in x and obtain 

T xt + v VT X + Vx -VT = AT X + f{T)T x . 



(4.13) 



Then, as in the proof of Theorem [| provided that A < C\ is sufficiently small, if we multiply 
( |4.13| ) by T x , and integrate, we obtain, using the Poincare inequality for T x : 



This implies that 



and hence 



lj t (\\TA 2 ) + C\\VT x f<\\vMVT x \l 



l|r,(.,t)|| 2 < ||t x (-,o)||V c '*/ a2 + c" / 

We integrate ( |4.14|) and obtain 



s)\\ 2 e' c '^ x2 ds. 



(4.14) 



t PS 



\\T x (;s)\\ 2 ds<C\ 2 \\T x (.,0)\\ 2 + C 




\u -,si e 



2~-C'(s- Sl )/X* dsids 



JO 



<C\ 2 \\T x (;0)\\ 2 + C\ 2 / M; S )\\ 2 ds. 



We use now the bound (|4.12| ) to get 

/ \\T x (-,s)\\ 2 ds<C\ 2 \\T x (-M 2 + —\H-M 2 + Cp 2 \ e \\T x (s)\\lds. (4.15) 
Jo a Jo 

Therefore there exists a constant Co > so that if ( |4.9|) holds, then 

\\T x {; S )fds <c\ 2 \\t x {-m\ 2 + — IK,o)|| 2 . 

a 

Then ((O]) and ( ^T2|) imply that 



H;s)\\ 2 ds < Cp 2 X 4 



r\ 2 

CA 2 ||T x (-,0)|| 2 + — 1|<,0)|| 2 



a 



a 



(4.16) 



C, 



<ci|t;(-,o)|| 2 + -|K-,o)|| 2 . 

O" 

The bounds ( ^4 . 1 4J ) and ( [4.16[ ) together imply in an elementary way that 

\\T x (-,t)\\ ^Oast -++oo. 
We multiply the vorticity equation ( |2.16| ) by uj and integrate in space: 



(4.17) 
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Then we have 



\co(;t)r<\\co(;r)r + 2ap I M;s)\\\\T x (.,s)\\ds 

1/2 / 



< \U;rW + 2ap 



\u>(-, s) II ds 



\T x (;s)\\ 2 ds 



1/2 



The uniform bound ([4.16|) implies that there exists a sequence of times — > +00 so that 
W^i'i T k) II — * 0. Then we have for t > r k : 



\\cu(;t)\\ 2 < |H-,r fc )|| 2 + 2ap 



s) \\ 2 ds 



1/2 



\T X (; s)\\ 2 ds 



1/2 



and hence ||o;(-,i)|| — > as t —>■ +00 because of ( |4.15| ) and ( 4.16 ). In order to prove ( |4.4| ) we 
note that ( [4.1 7|) implies that 



ta 



t-2 



|Vo;(-, S )|| 2 ^ < / ( Tp||cu(-, S )||||T,.(-, S )||^ + -||^(-,t 1 )|| 2 



't! Jtl * 

so that ( ^B| ) implies that for any e > we have for ti, ti > t\ + 1 sufficiently large we have 

\ \\Vuj(- } s)\\ 2 ds < e(t 2 -h). 
Jti 

Therefore, ( ft.llQ , (P-21 ) and the Cauchy-Schwartz inequality imply that 

V{t)<2 + - f \\w{s)\\ 00 ds + o{l)<2 , _ 

t Jo Vt \Jo 

and (|4.4|) follows. This finishes the proof of Theorem 



n / ft \ 1/2 

' ' "Vuj(-,s)\\ 2 ds ) +o(l) < 2 + e + o(l) 



5 Linear instability 

In the previous section we established the stability of planar fronts with respect to short 
wavelength perturbations. We analyze now the linear instability of planar front with respect 
to long wavelength perturbations. 

We perform a Galilean transformation z 1— ► z — v^t in ( |2.1| )-( |iO| ) following the flat front. 
We write T(x,z,t) = r{z — v^t) + 6{x,z — vot,t), and, with a slight abuse of notation, 
v(x,z,t) = v(x,z — vot,t). We also linearize equations (|2.1| )- (|2~^ ), dropping the terms that 
are quadratic in 9 and v. We obtain the linearized system 



an 2 

_ _ Vodz _ kA6 - fj'{r{z))9 = -wt' 



du 
~dt 



— v d z u — uAu + Vp = gA9e z 



and for the vorticity we get the equation 



dt 



— v d z uj — vAlu = gAd x 9. 



(5.1) 
(5.2) 

(5.3) 
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Using the same reaction units ( J2.11[ )-( PTT5D as before, equations (|5.1|), ( |5.3| ) become 

89 89 



and 



()i -2--A9- f(P(z))6 = -wP'(z) 



duj duj 89 
at 8z ox 



The nonlinear equations, in the same units and frame of reference, are 

du „ du . 8T 

— — + v ■ vlj 2- aAuj = ap— 

at oz ox 



and 



dT dT 

- + v .VT-2--AT = f { T). 



(5.4) 
(5.5) 

(5.6) 
(5.7) 



We wish to show that flat profiles cannot remain linearly stable with respect to sufficiently 
long wave infinitesimal perturbations. We start for simplicity with the case of the infinite 
Prandtl number. We retain (6.41) but we use the infinite Prandtl number limit of (157 



89 



which, together with (|2.18|) allows us to write the active scalar rule 

w = -p(d x ) 2 (-A)- 2 9. 
We express 9(x, z, t) in terms of its Fourier series: 

9(x,z,t)= J2 9k(z,t)e lkx . 



(5. 



(5.9) 



(5.10) 



The linearized temperature equation ( B.4[ ) transforms into 

^ - 2^ + (k 2 - d zz ) g k - f'(P)g k = pQKg k 



dz 



with k = ±y, ±2^, . . ., the operator K defined by the Fourier transform of 



Kg = k 2 (k 2 -d zz y 2 g 



and 



Q(z) = -P'(z) > 0. 
We will choose the smallest positive wave number from now on: 

2tt 



k 



A 



(5.11) 

(5.12) 
(5.13) 

(5.14) 



The operator K defined by ( |5.12 ) is given explicitly by a convolution with a positive function 



(Kg)(z) 



4k 



:i + k\z-(\)e- k ^g(C)dC. 
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(5.15) 



The expression ( |5.15[ ) is obtained by an elementary calculation, starting from 

i r°° 

{k 2 - g{z) = -J e- k ^g(Od( (5.16) 

and iterating. It is well known that the profile P is decreasing in the case of KPP, bistable 
and ignition nonlinearities |j24 |, so that the function Q is positive. Moreover 

Q{z) > ae- blzl (5.17) 

holds for all z, with a > and b > absolute numbers that depend only on the nonlinearity 
f(T). The positivity of Q and of the kernel of the operator K in the right hand side of ( [5.1 1 ) 
imply that the solution of the initial value problem ( |5.11| ) remains nonnegative if the initial 
data is nonnegative. Let us consider a function <p(z) which has the properties 

e~ k{zl < <f>(z) < Ce- klzl (5.18) 

with C > 1 and 

\(f)'(z)\ < Cke~ klzl , \<j>"(z)\ < Ck 2 e~ klzl . (5.19) 

Such a function is obtained for instance by gluing smoothly 3e l ~ z on z > 1 to 3e 1+z on 
z < —1 by a function bounded below by 2 on the interval [—1,1]. Then one rescales z i— > kz. 
We multiply the equation ( |5.11|) by <f>(z) and integrate. Using the properties of (j> an d the 
positivity of g k (z,t) we obtain 



d_ 

It 



J <f>(z)g k (z,t)dz > -M J <p{z)g k {z)dz + p J <P{z)Q{z){K g k ){z)dz (5.20) 



This follows from integration by parts, and bounds ( |5.19 ) on the first two derivatives of <fi by 
constant multiples of 0. The constant M < 2C(1 + k + k 2 ), so it is uniformly bounded for 
bounded k, for instance for < k < 1. Now we bound from below the term 

p J <f>(z)Q(z)(Kg k )(z)dz >^J J ae- b ^<P(z)e- k \ z \- k ^g k ((,t)d(dz (5.21) 

We used ( |5.17|) , (|5.15|) and the positivity of g k . We neglected the nonnegative term contributed 
by k\z — (\ in ( |5.15|) . Now we use ( |5.18| ): 



P 



J <P{z)Q{z){Kg k ){z)dz > ^ (| e-^ 2k ^dz^j (J </>(C)<? fc (C, • (5-22) 



We obtain the ordinary differential inequality 

i /*<*)»w* >- {mdm ~ M ) I * [z)3 * (z)iz - (5 - 23) 

and thus H^Hl^r) grows exponentially in time. Therefore we have the following theorem for 
the infinite Prandtl number case: 
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Theorem 4 Let P(z — 2t), u = be a planar, x -independent traveling front solution of the 
infinite Prandtl number Boussinesq system 

^ + v .VT-AT = f(T) (5.24) 

with 

- Av + Vp = pTe z , V • u = 0, (5.25) 

with front boundary conditions for T at z = ±oo, vanishing velocity at z = ±oo and periodic 
boundary conditions in x of period A. There exists a positive constant (3 > such that, if 

p\ > p, (5.26) 

then the solution P is linearly unstable. This means that there exist infinitesimal perturbations 
which grow exponentially, when viewed in a Galilean frame of reference moving with the 
traveling front. Their exponential growth rate is proportional to pX. 

A similar instability analysis may be applied to a convective system with the infinite Prandtl 
number equation ( |5.8|) replaced by the Darcy law that was studied numerically in (lCj . 



We return to the system ( |5.4j , |5.5| ) at a finite Prandtl number. We use the Fourier expansion 

w{x, z,t) = J2 w k(z, t)e 2mnx/x (5.27) 

nGZ 

and using ([2.18 ) we have, with k = 2irn/X, 



M*> *) = ~2 J e~ k]z ~ A Uk{z\ t)dz' (5.28) 

with the obvious notation for the Fourier coefficients of u. An elementary calculation involving 
solving the linear equation ( |5.5|) starting with the zero initial data u(x, z, 0) = 0, and using 
0, gives 

w k (z,t) = —j j g k (z 2 ,s)I k (z - z 2 + 2(t - s),t - s)dz 2 ds (5.29) 

where g k (z,t) are the Fourier coefficients of 9 as in ( |5.10|) and 

4( C , r ) = / e -*IC4W5^ c -l# ( 5. 30 ) 

J V 27T 

We now bound I k at this point as follows: 

4(C, t) > e- CTfcV / e-fclCI-fcMV^e-^-^L > 2 e - fc I^Erf(fcv / 2^) (5.31) 
J v 2n 



with 

Erf(a)= / ( 5 - 32 ) 



2ti 
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For a > 1 we will use just a piece of this integral, 



r 2a „a du 1 

Erf (a) > / e~~^= > —^=t 



'2tx V27T 



and for a < 1 we write 



Erf (a) > Erf(l). 



This allows us to deduce, for later use 



4(C,r)>2(l 



f 2ir' 



>-fc|Cl e -2 CT fc 2 7 



(5.33) 
(5.34) 

(5.35) 



The evolution equation for follows from (5~4") 



lit ~ 2 ~iz + ^ 2 ~~ dzz ^ 9k ~ f'( P )9k = (TpQw k 



(5.36) 



with Wk computed using Q5.29|) and with Q = —P' as before. Because of the explicit expression 
( 5.29|) we see that if the initial data gk(z, 0) is real and positive, then Qk(z, t) remains real and 
positive. We take such initial data. Then, using ( |5.35|) and the triangle inequality we get, 



Wk(z,t) > (1 -=)pake 

v2ir 



-k\z\ I e -2(t-s)k{ak+2) 



gk{z 2 ,s)e 



- klz2l dz, 



ds. 



(5.37) 



As in the case of the infinite Prandtl number, we take a function <j>{z) with properties ( |5.18| ), 
( |5.19| ), multiply ( |5.36|) by and integrate. Let us denote 



Y(t) = / 4>{z)g h {z,t)dz. 



(5.38) 

Integrating by parts we get 

dY r l 

—- + C 3 (l + k + k 2 )Y > C^apk \ e- 2{t - s)k{(Tk+2) Y(s)ds. (5.39) 
dt J 

with C3 and C4 positive numbers that depend on the properties of <fi. Consider also 

Z(t) = e 2k{ak+2)t Y(t). (5.40) 
Multiplying by e 2fc(CTfc+2) *, we get from ( |5.39| ) 

+ (3Z > a [ Z(s)ds (5.41) 
Jo 



dZ_ 
~dt 



with 
and 



(3 = C 3 (l + k + k 2 ) - 2k{ak + 2) 
a = C^apk. 



(5.42) 
(5.43) 
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It is clear that solutions of the differential inequality (|5.41 ) are larger than solutions of the 
ODE 

dy 



ijt +(3y = a I y(s)ds (5.44) 



o 



with the same initial data. This ODE is solved differentiating one more time, 

— + p- - ay = (5.45) 

and seeking solutions of the form e M *. Equation ( |5.45|) always has at least one exponentially 
growing solution because a > (irrespective of the sign of 0). The general solution of ( |5.45| ) 
is 

y(t) = y 1 e t ^ +y 2 e t »- (5.46) 

with 

fi± = \ (~P ± v / /5 2 + 4a) . (5.47) 

The solutions of ( |5.44[ ) correspond to the linear subspace for which ^(0)+/3y(0) = 0. Choosing 
without loss of generality the coefficient of the growing exponential y\ = 1, we deduce the 
relation 

p + + (3 + y 2 (p.+(3)=0 (5.48) 

In order to have exponentially growing positive solutions we need the initial datum for ( |5.44 ) 
to be positive. This initial datum is 1 + y 2 and is computed using ( |5.47| ) and (|5.48|) 



l + y 2 = , 5.49 

so it is always positive. Starting with this initial condition (or any positive multiple thereof) 
we get exponential growth for Z(t). This will imply exponential growth for Y(t) if 

H+ > 2k(ak + 2). (5.50) 

This turns out to be the condition 

C i ap>2C 3 (l + k + k 2 )(ak + 2) (5.51) 

or, putting C = 

p> + (l + k + k 2 ) (5.52) 



Theorem 5 Planar reactive Boussinesq fronts ( \5.(\ , |5. 7 ) are linearly unstable to large wave- 



length perturbations whenever the local Rayleigh number p based on the laminar front thickness 
is large compared to the inverse of the Prandtl number, 

2C 

P> — • 
a 

Perturbations with wave numbers k satisfying ( \5. 53{ ) grow exponentially in a frame of reference 
moving with the planar front. The growth rate is proportional to \fo~pk. 
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Remark. The exponential growth rate proportional to the square root of the wave number is 
a signature of the Rayleigh- Taylor instability, operating here only at large scales. When the 
wave length of the initial perturbation is decreased to a length comparable to the thickness 
of the planar front, the perturbation decays in time. 

6 Conclusions 

We proved that the reactive Boussinesq system in a strip has front solutions with bounded 
speeds, if the nonlinearity is of concave KPP type. The bounds grow with the aspect ratio, 
but nevertheless, front acceleration does not occur in this system. The concavity of the 
nonlinearity was used in the proof. We do not know at present whether the divergence of the 
bound with aspect ratio is sharp. For small aspect ratios and for small Rayleigh numbers, 
the only traveling modes are planar, and all front-like solutions become planar. This result is 
expected to be valid in more generality than proved here. For large enough Rayleigh numbers, 
if the aspect ratio is large, then the planar fronts lose stability to long wave perturbations. 
The instability is of Rayleigh- Taylor type. This again is expected to hold in greater generality 
than proved here. The results proved here agree qualitatively with the recent numerical study 

H- 
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